Cosmological parameter dependence in local string theories of structure formation 
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We perform a new and accurate study of the dependence on cosmological parameters of structure 
formation with local cosmic strings. The crucial new ingredients are the inclusion of the effects 
of gravitational backreaction on the evolution of the network, and the accurate evolution of the 
network through the radiation to matter transition. Our work re-iterates the fact that expanding 
Universe numerical simulations only probe a transient regime, and we incorporate our results into 
the unequal time correlators recently measured. We then compute the CMB and CDM fluctuations' 
power spectra for various values of the Hubble constant Ho and baryon fraction SI5. We find that, 
whereas the dependence on Qb is negligible, there is still a strong dependence on Ho. 
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I. INTRODUCTION 



Modern theories of structure formation are based on two paradigms: cosmological inflation [Q and topological defects 
||. Inflationary scenarios produce structure from vacuum quantum fluctuations, and have been well studied. Hence 
upcoming high precision experiments aimed at measuring the Cosmic Microwave Background (CMB) anisotropy and 
£f~) • Large Scale Structure (LSS) power spectra should also provide information on the model's free parameters. Parameters 
internal to the theory are, for example, the primordial tilt, the perturbations' amplitude, or the percentage of tensor 
ON ■ modes. Others describe the current state of the Universe, and include the ratios to the critical density, of the matter 
' density, f2, of the vacuum energy, of the baryon density, fit,, and the Hubble constant Hq. 

Progress in defect theories has been slower. Here, as the Universe cools down, high temperature symmetries are 
Cln' spontaneously broken, and remnants of the unbroken phase may survive the transition p|. These "vacuum" defects 
later seed fluctuations in the CMB and LSS. The defect evolution is highly non linear, thereby complicating the 
computation of these fluctuations. However, recent breakthroughs in method and computer technology have allowed 
unprecedented progress although the results for a standard combination of parameters (fi = 1, fl\ = 0, 

fib = 0.05, Hq = 50 km/sec/Mpc) were discouraging. Tuning parameters to fit the data has not been thoroughly 
attempted, although initial indications are that it is possible to improve the fits considerably Also, combinations 
• i-h . of strings and inflation with the standard parameters seem to fit the data much better than each separate component 

Unlike inflation, defect theories only have one internal parameter: the symmetry breaking energy scale (which is 
also the string mass per unit length /i for local cosmic strings). Nonetheless, quantities like O, O & , S1a, or Hq are also 
free parameters in defect theories, although it has been suggested in the context of global theories that the dependence 
on these is much weaker ||. The purpose of this paper is to reexamine this statement, in the context of local cosmic 
strings. 

Before such a study can be attempted a number of aspects of the work on strings jl[j6) need to be refined. The 
effect of the radiation to matter (RM) transition, and the transition to curvature or vacuum domination upon the 
network need to be established, if sensitivity to the cosmological parameters is to be fully considered. However, these 
effects were neglected in the recent highest accuracy calculations of structure formation with local cosmic strings [[| , 
which make use of unequal time correlators (UETC). Here we will use the 3-scale model of Austin, Copeland and 
Kibble flT^-|l4|] to evaluate those effects in a flat Universe with no cosmological constant. A further key ingredient of 
this model is the inclusion of gravitational backreaction (GBR). These results are incorporated into the UETCs of ||, 
and the CMB and LSS power spectra are obtained on a grid of values for fi^ and Ho. We argue that this combination 
of techniques provides the most accurate calculation of these spectra allowed by current computer technology. 

This paper is set up as follows. First, in section [nj we discuss the aspects of the 3-scale model which are relevant 

[|. The resulting 
Conclusions are 



to this work. In section III we then explain how these results are incorporated into the UETCs of 



CMB and LSS power spectra, and their dependence on Ob and Hq are then presented in section |V 
given in section 
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II. THE 3-SCALE MODEL 



The 3-scale model |T^,[T3|] provides an analytic description of the string network in terms of three physical lengths 
£, £ and £. The mean string separation is £ so that the energy density in strings is p = fj,/^ 2 . (Renormalisation 
effects are encoded in £ and not /x.) The mean correlation length along the strings is £(7^ £). Finally the small scale 
structure is described by £ which is effectively the interkink distance. These scales are implicitly functions of each 
other, of time t and of parameters describing the effects of expansion, gravitational radiation, GBR, loop formation 
and intercommutation. They also determine the root mean square string velocity v rms . In this paper we have solved 
the evolution equations for d^/da, d^/da and dC,/da in terms of the scale factor a and the variables 7 = 1/H£, 
7 = 1/H£ and e = 1/HC, where H is the Hubble parameter. We do not write out these rather long equations here: 
they follow directly from those given in JL2J by a simple change of variables from t — > a. Of particular interest is the 
effect of GBR, the behaviour of the network across the RM transition, and the question of whether or not the system 
has reached a scaling regime (e, 7 and 7 all constant) today. 

The 3-scale model could be criticised for containing many undetermined parameters. However, some of these — such 
as the infinite string intercommutation probability Xi the chopping efficiency c which gives the rate of loop formation, 
and TG/i the rate of gravitational radiation, also appear in 1-scale models, for example ]l5| ]. Other parameters are 
tightly constrained by numerical simulations of cosmic string evolution (see below). Indeed, the only undetermined 
parameters are k, which describes the excess small-scale kinkiness on loops compared to long strings, and C which 
determines the rate at which GBR smooths the small-scale kinkiness; d(/dt\ GBR = CTG/i. Neither quantity is to 
be found in current numerical simulations nor in other analytic models. One further point; we assume here that 
gravitational radiation and not particle production is the dominant source of energy loss [ |16| . This implies that 
C > Ccrit > k for C cr it a given constant of order 1 |L2| . The reader is referred to jl^] for a recent discussion of the 
effect of different decay mechanisms on the observational consequences of strings. 

Some important features of the solutions for 7, 7, e and v rm s are shown in Fig. [l] for different values of C. One of 
these is that there are two different scaling solutions. The first, 'transient' solution, has constant values of 7t r , 7t r 
and u'^ {independent of C), but e is growing: small scale structure is building up on the network. We believe that 
this corresponds to the scaling solutions of numerical simulations which have no GBR (so (7 = 0). From the evolution 
equations one can obtain expressions for 7 tr , jf r and i>*„ s [ fT^ and thus bound them by data from simulations, which 
in turn constrains many of the 3-scale model parameters as claimed above. We comment that the duration of the 
transient regime does depend on C, ending when e ~ 1/TGfiC. Thus if, as in numerical simulations, we were to 
set C — it would last indefinately, with e — > 00 eventually causing the strings to disappear (7,7 —> 0). Finally, 
recall that in Ref Jl6| particle production was found to be the dominant source of energy loss from the network. 
The analogue here correpsonds to having k > C cr it for which the evolution equations lead to a scaling regime with 
7 ~ e ~ 7; a 1-scale model |Q. 

Fig. ^ shows that following the transient regime is a new 'full' scaling regime in which all the variables reach new 
constant values, apart from a period between — 5< ln(a/a eq )<10 (including recombination and even today) resulting 
from the RM transition. It appears that the string network might only just, if at all, be reaching a scaling solution 
today. Indeed, across the transition, different scales seem to reach scaling at different rates which are also dependent 
on C. We also see that as C increases £ increases since kinks are smoothed out more effectively; v rms is also increased. 
It is this full scaling solution, and not the transient one that is clearly the relevant one for structure formation. 
However, we noted above that C is so far unknown. We have therefore chosen to work in the limit where the scaling 
solutions are as insensitive as possible to C. From Fig. |l| this corresponds to large values of C for which the decrease 
in physical energy density across the radiaton-matter transition is also largest. In that limit, the solutions plotted in 
Fig. |l] are also insensitive to the precise value of k(< C cr u) — the exception is e which decreases as 'particle production' 
(k) increases. 

III. METHOD OF DETERMINATION OF THE POWERSPECTRA 

We have used the full scaling solutions including GBR with large C to update the results of ||, where the power 
spectra in string induced CMB fluctuations (Cg) and LSS (P(k)) were computed following methods first proposed 
in jj). Of key importance is to obtain from simulations the 2-point functions of the defect stress-energy tensor Q^. 
These so-called unequal time correlators (UETC) are defined as (0 /Jt/ (k, r)0* /3 (k, t')) = C fJil/ _ a fj(k, r, r') where k is the 
wavevector and r and r' are any two (conformal) times. UETCs contain all the information required for computing the 
power spectra in defect induced fluctuations. They fall off sharply away from the equal-time diagonal (r = r'), and are 
tightly constrained by requirements of self-similarity (or scaling) and causality. These properties allow extrapolation 
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FIG. 1. Evolution of 7, 7, e and v rms for different values of C and Ho = 50 km/sec/Mpc. Other parameters are k — 0.05, 
X = 0.1, c = 0.15 with TGfj, = 10~ 8 to emphasise the transient scaling regime (see text). 




FIG. 2. Plot of a, the comoving energy density in strings multiplied by the conformal time squared 
(a — Np c T 2 — NjiH 2 ^^^ 2 where the normalisation TV is such that a = 1 deep in the radiation era). Different curves 
correspond to different values of C and other parameters are as in the previous caption. 
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FIG. 3. The CDM power spectrum in sCDM scenarios (top curve), and local cosmic strings with large C. We consider the 
result of a Boltzmann code calculation, with and without the effects upon the network of the RM transition included, and also 
the results of a Albrecht and Stebbins approximation. The latter underestimates the normalization and distorts the shape of 
P(k). 



of measurements made over the limited dynamical range of simulations, thus giving the defect history over the whole 
life of the Universe. 

However, it is clear that scaling breaks down near the RM transition (Fig. |l|). This may be incorporated into the 
UETCs by measuring radiation and matter epoch UETCs, transition UETCs, and then interpolating ||. Such a 
procedure was applied to global defects, but cannot be followed with the simulations of j^], in which the expansion 
damping effects upon the network are neglected. Nonetheless these simulations are still the only available complete 
source of local string UETCs (see however ) . The 3 scale model of jl2j , while allowing the most rigorous study 
to date of departures from scaling, does not provide sufficient information for constructing UETCs. 

We therefore decided to adopt a hybrid approach, and incorporate the information supplied by the 3-scale model of 
M| into the UETCs of || ■ The UETCs are renormalised with a non-scaling factor describing the change in comoving 
string density across the RM transition, as determined from the 3-scale model. That is we make the replacement 

C/wrfi^y) ->■ C-^u,a0(k,T,T')a(T)a(T') (1) 

where a(r) measures the change in the defect scaled comoving density. More specifically, a(r) = Nt 2 p c = Nfi(HaTj) 2 
where N is a normalisation factor chosen such that a(r) — 1 deep in the radiation era, and p c is the comoving energy 
density in srings. Note that a depends on 7 and hence (see figure JO) on C: this C dependence of a is plotted in 
figure ^. The procedure given in equation (|I|) was also followed in [j7|, in the context of a semi-analytical model for 
string induced fluctuations. The rationale lies on the fact that UETCs are proportional to the square of the defect 
comoving density, if one assumes no correlation between individual defects (as in the model of Q). The presence 
of inter-defect correlations, non-negligible for small wavelengths, means that the shape of the correlators could in 
principle also change with 7. However, it is difficult to extract this information from the 3-scale model, and we expect 
the effect upon the CMB and LSS to be small. Naturally all UETCs (scalar, vector, and tensor) are multiplied by 
the same time-dependent factor. This then propagates to the various eigenmodes. 

It is crucial to notice that it is the comoving density, and not the physical density that enters the above scaling 
argument. The conversion factor between physical and comoving densities changes appreciably from the radiation to 
the matter epoch. Indeed for large C it accounts for most of the drop in defect physical density, and as a result the 
strings' comoving density (or a(r) plotted in figure |J) drops only by about 20% for large C. For small C the comoving 
energy density actually increases as can be seen from figure |^. 

We reran the codes used in H with these modifications (working throughout with large C), and with two purposes 
in mind: to compare the effects of neglecting the RM transition with approximations used by other groups, and to 
survey parameter space in cosmic string models. 
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FIG. 4. The CMB power spectra predicted by cosmic strings for different values of the Hubble constant. We have plotted 
(£(£+ l)Ce/2n) 1 / 2 in fiK, and superposed several experimental points. The effect of varying Qb is negligible, and so for clarity 
in this picture it was kept fixed. 



IV. RESULTS 



We found that the effect of neglecting the RM transition is not very large, and certainly much smaller than the errors 
induced by using the Albrecht and Stebbins (AS) approximation pl| (or a variation thereof As suggested in 

||, Gfi normalized to COBE increases by a factor of the same order as the decrease in string density in the matter 
epoch. With standard parameters (CI = 1, f^A = 0, f2& = 0.05 and Hq = 50 km/s/Mpc) we found that G[i increases 
from 1 x 10 -6 to 1.2 x 10 -6 The COBE normalized CV spectrum, in this cosmology, hardly changes at all. The 
bias in spheres with radius 100 Mpc hr 1 changes from 4.9 to 4.2: the inclusion of the RM transition improves the bias 
problem of cosmic strings but not by much. The effect on the shape and normalization of P(k) is plotted in Fig. |^. 
By contrast, using the AS approximation instead of a Boltzmann code (or one of its truncations [fl^lf) induces large 
uncertainties in normalization and shape of P(k). (We obtained the AS curves of Fig. ^| by replacing the variables £ 
and x °f 01 with the £ and £ of the 3-scale model respectively.) 

Although a small effect, only once the RM transition has been properly taken into account can we start exploring 
dependence on cosmological parameters. We considered fib and h, and solved for the Ci and P(k) on a grid with 
0.4 < h < 0.8 and 0.01 < flbh 2 < 0.015. The latter satisfies nucleosynthesis constraints, the former accommodates 
currently accepted ranges of the Hubble parameter. In Figs. |i| and || we plot Cg and P(k) extracted from this grid 
to illustrate the observed variations. The Cg spectrum of anisotropics induced by strings seems insensitive to fib, but 
depends sensitively on h (Fig. |J). In particular, local cosmic strings, unlike global defects, exhibit a Doppler peak, 
but only if the Hubble constant is not too large. The results of j5j are therefore not generic, but depend on the choice 
of h; indeed they are consistent with the results of j2^] in which strings are directly coupled to a Boltzmann code. In 
that work it was found that strings do not have a Doppler peak, with h = 0.8. Notice how the dependence of the Ce 
on h is non-linear, with a much larger jump from h — 0.7 to h — 0.8 than, say, from h = 0.6 to h = 0.7. The effect of 
both h and Qb upon the CDM power spectrum is very small if P(k) is plotted using units of /i/Mpc for wavenumbers 
(Fig. ||). This is in agreement with small k analytical predictions (eg. p3[), but breaks down at large k as can be 
seen by the diverging curves in Fig. ||. The string bias problem at 100 Mpc h^ 1 is clearly not solved by changing h 
or fib (but see [^|-|| for other solutions). Finally, recall that we have worked with large C. For small C"s the strings' 
comoving energy density increases across the RM transition and the bias problem is in fact worse. 



V. CONCLUSION 



To summarize, in this paper we have studied the impact of GBR on cosmic string network evolution, and properly 
accounted for its properties during the RM transition. Not only did we suggest that this departure from scaling lasts 
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FIG. 5. The effect of changing h and fl^h 2 upon the power spectrum P(k) in CDM. 

longer than that obtained in other models, but we re-iterated the existence of two different scaling solutions. The full 
scaling solution is the relevant one for structure formation and cannot be observed in numerical simulations. 

These results, with large C, were then used to update the calculations of ||, and cut a first, safe slice through 
parameter space in string scenarios. Dependence on h and fib was considered, the conclusions being that the CMB 
results are sensitive to h but less so to fib- If forthcoming experiments were to reveal a Ci spectrum without 
secondary Doppler peaks — a tell-tale signal of a defect theory Q — then we would lose ability to determine fib, since 
information on fib is hidden in the secondary Doppler peaks which are generically erased by incoherence, in defect 
scenarios. However, for h, our conclusions are more positive. Local strings exhibit a Doppler peak whose height we 
have demonstrated decreases quite sharply with an increasing Hubble constant. 

Larger regions of parameter space require surveying. In particular as pointed out in refs strings with a 

cosmological constant are a promising combination. However we feel that all treatments, including our own, become 
too qualitative to be reliable if curvature or A are present. Considerable analytical and numerical work with fl^ ^ 0, 
and in open models, is required before a more quantitative discussion is possible. The first step in the context of the 
method we are proposing here would be to include A into the 3-scale model. 
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